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1. Introduction 

All neutral mesons — K Q , B°, B s , D° — have been observed to oscillate from particle to antipar- 
ticle. The oscillation frequency AM tests the Standard Model's pattern of flavor violation. The 
phenomenology is especially simple for neutral B mesons (normal and strange), because the flavor- 
changing dynamics play out predominantly at distances much shorter than the scale of QCD. In the 
case of the B mesons, the width difference Ar of the two propagating eigenstates also arises pre- 
dominantly at short distances. It is especially intriguing (at least for now), because measurements 
of AT S and the CP phase <j> s of the B s are in imperfect agreement with the Standard Model |[l], Q|. 

Neutral B mixing stems from AB = 2 flavor-changing transitions. In the Standard Model 
these arise first at the one-loop level, so non-Standard contributions are conceivably of com- 
parable size. The observables are then (approximately) AM = 2\M\2\, AT = 2|ri2|cos0, and 
= arg (— M12/T12), where M12 and Ti2 are the off-diagonal elements of the mass and width 
matrices of the two-state systems: 

M12 = ^^(y^) 2 5o(m f 2 /Mi,)j] fo (Ai)(B|fe7^^|B} +BSM, (1.1) 

r i2 = -Prr- [G(V,ll)(B\q L Y^bq L fb\B)+G s (V,ll)(B\q L bq L b\B}] +BSM, (1.2) 

where V is the CKM matrix, and Sq, T]^, G, and G$ are short-distance effects, computed in elec- 
troweak and QCD perturbation theory. Contributions beyond the Standard Model ("BSM") are 
not written out explicitly. Because of the V — A structure of the electroweak interaction, only the 
left-handed (light) quark field q\ = qh(l + 75) appears. 

The remainder of this paper is organized as follows. Section |2] constructs lattice operators 
with staggered light quarks and Fermilab heavy quarks, corresponding to the 4-quark operators in 
Eqs. ( |P| ) and (L2). (The construction suffices for any light quark with chiral symmetry and heavy 



quark with heavy-quark symmetry.) We give a status report of our numerical results in Sec. ||. 
Section |] summarizes and presents some of our plans for the future. 



2. Short-Distance Matching 

To compute the hadronic matrix elements in Eqs. Ql.l[ ) and ( |1.2| ), one has to derive an expres- 
sion in lattice gauge theory that approximates well qLj^bqij^b and qibqib. The lattice operators 
can then be computed, and the numerical and other uncertainties estimated, to determine M\i and 
Ti2. Similar operators appear BSM, for which the following derivation serves as a template. 

For the light valence quark we take naive asqtad propagators 

(Y(x)Y(y)) u =a(x)Q.- l (y)(x(x)x(y))u, (2.1) 

where % is the one-component staggered fermion field; T is a 4-component naive field, and (• • -)u 
denotes the fermion average in a fixed gauge field U. For the heavy quark we use 

¥=[l+d 1 (m a)Y-D]\if, (2.2) 

where y is the fermion field appearing in the Fermilab action Q or an improved action with the 
same design features [[J. 
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We aim to construct lattice operators Q and Q$ such that 

Q = qLY^bq L Y^b + 0(a 2 ), (2.3) 
Q s = q L bq L b + 0(a 2 ), (2.4) 

where = means "has the same matrix elements as." Here the 0{a 2 ) term depends on mba. As long 
as one retains small corrections to heavy-quark symmetry, it remains bounded even as n%a — > °°; 
as long as certain Dirac off-diagonal improvements are consistently introduced [Q, ffl], they vanish 
as a — ► 0. These two elements are the essence of the Fermilab method. 

Our construction starts with the lattice operators TlJ^ x ¥YlY^ x ¥ and Ti x ¥Tl x ¥. According to 
the HQET theory of cutoff effects []5], [g, 0], these lattice operators can be described by 

TlY^l^ = 2C M q L J li hMq L th^ + 2SC M q L h^q L h^ + £sj at ^ + • • ■ 
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(2.5) 
(2.6) 



i=\ 



where h^' are the heavy-quark fields of the heavy-quark effective theory (HQET), satisfying 
h& = i(l±7 4 )/j( ± ). The sums are over five dimension-7, AB = 2, four-quark operators, similar 
to those written out, but with an extra derivative. The series continues with operators of dimen- 
sion 8 and higher. On the right-hand side of Eqs. ( |2~31 ) and (^1]) the operators are to be understood 
with some continuum regulator and renormalization scheme. Discretization effects are lumped into 
the short-distance coefficients , 5C|^, and which depend on the couplings of the lattice 
action, as well as the lattice spacing a and the (renormalized) gauge coupling and quark masses. 

The next step is to note that the target operators have a completely parallel description in 
HQET, namely 

5 

qLYnbq L fb = 2Cq L y^ qrfh^ + 28Cq L h^q L h^ + + • • • , (2.7) 



q L bq L b = 2dC s q L Y fl h (+ k L Y^h^ + 2C s q L h^q L h^ + + 



i=\ 



(2.8) 



where the (continuum HQET) operators on the right-hand sides of Eqs. (|2.7j) and (2.8) are precisely 



the same as those on the right-hand sides of Eqs. (2.5) and (2.6). The coefficients differ, however, 



because the lattice does not appear on the left-hand side of Eqs. (2.7) and (2.8) 



With Eqs. ( |2.5[ )-(]2.8[) the desired construction of Q and Q$ is immediate: 

q s = Zst L wr L ^+sz s f L r tl wr L f^+Y i b Si Q i , 



(2.9) 
(2.10) 



where the <2r are lattice discretizations of the =2,-, such that Qj 
algebra then shows that if 



C\f&j + dimension 8. Simple 



Z 
8Z 



[ccf - scdcf] I [c lat c 

[5C-Z5C lat ] /cf, 



lat 

S 



8C l3X 8cf] 



[Bj-ZBf- 



5ZBfj\ C lat 



-l 

ji ' 



(2-11) 
(2-12) 

(2-13) 
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then Eq. (2.3) is satisfied. Similar expressions exist for Z$, 8Z$, and bsi, such that Eq. (2.4) is 



satisfied. From the structure of Eqs. ( |2.11| )-( |2.13| ) it is clear that the regulator and renormalization 



scheme dependence of the HQET drops out of Zm, 8Zi s \, and b^. 

Let us close this section with a few remarks. The enumeration of the operators and further 
operators of dimension 8, is an easy extension of Ref. [g]. In perturbation theory Crg\ (SC^ and 
the Bj) start at tree (one-loop) level, but they could also be determined nonperturbatively, adapting 
schemes such as that of Ref. [Jsj] . Because of the way Fermilab lattice actions are constructed [|], 
|J, starting with Wilson fermions, one has lim a ^oC lat = C, etc., without fine tuning. (In lattice 
NRQCD this is possible only with fine tuning.) Although our derivation hinges on the HQET 
description of cutoff effects, one could also (for m^a -C 1) use the Symanzik theory; the results for 
Z( 5 ), 5Z( S ), and b( S y would be the same. 

We have embarked on a one-loop calculation of Zrg\ and 8Zi s \. At present they are being 
checked. As with currents [§,01, it may prove prudent to write 

z (s) = z v bb Zv qq P(s), (2- 14 ) 

where Zy bb and Zy qq are nonperturbatively determined matching factors for the vector current. The 
remaining factor pi s \ could have a tamer perturbative expansion, because of cancellation among 
diagrams. We do not expect the cancellation to be as good as in the case of currents, because 
4-quark operators have new diagrams in which a gluon is exchanged from one bilinear to the other. 



With the rotation of Eq. ( |2.2[ ), the b^, in Eqs. ( J2.9[ ) and ( J2. 10| ) are of order a s and are not 
available. The calculations of the 4-quark operator matrix elements described below thus have 
discretization errors of the form 

»«*<*> ...a „ (2.15, 



{qLYiibqiy^b) 2(1 +m Q a) 

dim 8 ops ~ a 2 A 2 f(moa), (2.16) 

where the mass dependence of the B/ s ^ is an Ansatz with the correct asymptotic behavior as 
mo — > 00 and as m§a — > for the Fermilab action. The functions f(moa) multiplying the 0(a 2 ) 
discretization effects are known for the Fermilab action |9|]. 

3. Long-Distance Matrix Elements 

To compute the matrix elements we use a data-object called the open-meson propagator [[R]]. 
Valence quark propagators are started at an origin (xo,?o), where the 4-quark operator sits, out 
to all (x,t). Since, for this problem, we are interested only in zero-momentum pseudoscalars, at 
each t the Dirac indices are contracted with 75, and this contraction is summed over all x. On the 
other hand, M12 and Ti2 require two (several) Dirac structures in (beyond) the Standard Model. 
Therefore we leave the Dirac and color indices free at (jco,?o)> writing out one 12 x 12 x N4 data- 
object per configuration, where N4 is the total number of time slices. Three-point functions are 
formed by contracting open-meson propagators at times ?, and tf with the Dirac structure of each 
4-quark operator. Two-point functions from to to t are used to normalize the matrix elements and 



to provide a cross-check with our separate calculations of B-meson decay constants [11]. 
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Our calculations are carried out on several ensembles of lattice gauge fields with a realistic 



sea of 2+1 flavors, made available by the MILC Collaboration [|12j, |13p. The ensembles used here 
are listed in Table [j] together with the valence quark masses. The sea quarks are simulated with 
the asqtad action for staggered quarks, and with the fourth-root procedure to reduce the number of 
species from 4 to 1 . 

To discuss the analysis, it is helpful to introduce some notation. The four-quark matrix ele- 
ments are written 

(B°|T L 7^ l7 ^|B°) = lM Bc $, (3.1) 

where the quantity fi q is well-behaved in the heavy-quark limit. We extract fi s and /3 C / from 2- and 
3-point functions. With staggered valence quarks these correlators have contributions from wrong- 
parity states with time dependence (— We are careful to disentangle these states. To isolate 
the ground state we use Bayesian fits, varying the number of states. 

We then carry out a partially-quenched (i.e., m q and mi varying independently) chiral extrap- 
olation of j3 9 /j8 $ to obtain j3rf/j8 $ , using rooted staggered chiral perturbation theory for fi q JM) , p"5| ]. 
With more valence masses than sea masses, the effects of partial quenching constrain the param- 
eters of #PT more stringently than would unitary (m q = mi) data alone. Fitting the ratio j3^/j3 s 
yields smaller statistical errors than fitting r 3 J 2 f5 q directly. We also carry out a chiral extrapolation 

3 /2 

of r x fi s , which is mild, because it depends only on the sea masses (ami,amh). 

In the phenomenology of B-B mixing it is conventional to write the matrix element as 

(B q \q Lrix bq L fb\B q ) = lf 2 B M 2 B B Bq . (3.2) 



Neglecting Z — 1 and 8Z in Eq. ( j2T9[ ) one sees that p q = f Bq JM Bq B Bq . Of special importance is 

% = fBB^/fsXi 2 = (M B jM B f/ 2 (P s /M, (3-3) 
where, again, the right-most expression neglects Z — 1 and SZ. We use the experimentally measured 

^ 1/2 

meson masses and our chirally extrapolated j8 s and jS^ / j8. v to obtain f B B B ' and t, . The light-quark- 
mass dependence is shown in Fig. [j] Further plots can be found in Ref. [16]. 



A preliminary, but comprehensive, error budget is given in Table ||. The B*-B-n coupling g B * B7[ 
enters the expressions for the chiral extrapolation. The data are not precise enough to determine 



a (fm) 


Lattice 


N confs 


Sea (ami,amh) 


Valence am q 


0.12 


24 3 


x64 


529 


(0.005,0.05) 


0.005, 0.007, 0.01, 0.02, 0.03, 0.0415 


"coarse" 


20 3 


x64 


833 


(0.007,0.05) 


0.005, 0.007, 0.01, 0.02, 0.03, 0.0415 




20 3 


x64 


592 


(0.01,0.05) 


0.005, 0.007, 0.01, 0.02, 0.03, 0.0415 




20 3 


x64 


460 


(0.02,0.05) 


0.005, 0.007, 0.01, 0.02, 0.03, 0.0415 


0.09 

"fine" 


28 3 
28 3 


x96 
x96 


557 
534 


(0.0062,0.031) 
(0.0124,0.031) 


0.0031, 0.0044, 0.062, 0.0124, 0.0272, 0.031 
0.0031, 0.0042, 0.062, 0.0124, 0.0272, 0.031 



Table 1: Input parameters for the numerical calculations. The lattice spacings listed are approximate 
mnemonics. The heavier sea mass ;«/, is close to the strange mass, which then is subject to retuning a poste- 
riori, yielding the last value of am q for the coarse ensembles. 
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gB*B%, so it must be set with a prior distribution in the chiral fits. A range that encompasses 
phenomenological and quenched lattice estimates is gB'Bx = 0.35 ±0.14. The error in Table | 
corresponds to this range, while the prior width in the fits is ±0.28. 

1 12 

Until the perturbation theory has been checked, we prefer not to report a value for fB s Bg s . 
The matching corrections nearly cancel in the ratio p q /p s ; the results with and without Z — 1 and 
SZ are nearly the same, as shown in Fig. [I]b. With the error budget discussed above we find 

E, = 1.205 ± 0.037 stat ±0.034 syst , (3.4) 

unchanged since Lattice 2008 [|I"5|]. 



4. Future Prospects 

When the perturbative matching has been completely checked, we will be in a position to 
present final results. We can also compare different strategies, in particular, whether the perturba- 



tive expansion seems to work better for p/ s \ or Zt$\ (cf. Eq. ( |2. 14[ )). 

In the longer term, we plan to obtain results for 4-quark operators that enter beyond the Stan- 
da**^ ^Ar+AcA Ci . ft li .-.-w-i *~\ fUa r~* < ,-i i-w-i Ul m *"*VV * lw " 1 „4-„4-:„4-: — 4-1 4-1 — 



(a! 



Preliminary 



0.02 0.03 0.04 0.05 0.06 0.07 



1 II 

Figure 1: Light-quark-mass dependence of fB s B^ s and ^ . The curve in the right plot is a fit to all partially- 
quenched data, not just the shown unitary data. 



S 

pa' 



(b 



Preliminary 



# a-0.12fm.FNAL/MlLC 

V a=0.09 fm, FNAL/MILC 

+ a=0.12 fm, FNAL/MILC, renormalized 

▲ a=0.09 fm, FNAL/MILC, renormalized 

— chiral extrapolation (after a— and fix m ) 




0.04 
r ,fn,. ,. 



Source 


Av 


% 


5 


Statistics 


2.7 


4.0 


3.1 


Scale (n) 


3.0 


3.1 


0.2 


Sea and valence quark masses 


0.3 


0.5 


0.7 


&-quark hopping parameter 


<0.5 


<0.1 


<0.1 


^PT + light-quark discretization 


0.4 


2.5 


2.8 


gB*Bn 


0.3 


0.6 


0.3 


Heavy-quark discretization 


2 


2 


0.2 


Matching (perturbation theory) 


~ 4 


4 


<0.5 


Finite volume 


<0.5 


<0.5 


<0.1 


Total 


6.1 


7.3 


4.3 



Table 2: Preliminary error budget. Entries in percent. 



6 



B®-B® Mixing and Matching Andreas S. Kronfeld 



current project at a = 0. 12 and 0.09 fm, but also extend to smaller lattice spacings, a = 0.06 and 
0.045 fm. New runs with higher statistics and five lattice spacings (also 0.15 fm) are underway. 

References 

[1] A. Lenz and U. Nierste, Theoretical update ofB s -B s mixing, JHEP 0706 (2007) 072 
[arXiv:hep-ph/0612167]. 

[2] M. Bona et al. [UTfit Collaboration], First evidence of new physics in b <-» s transitions, 
arXiv:0803.0659 [hep-ph]. 

[3] A. X. El-Khadra, A. S. Kronfeld, and P. B. Mackenzie, Massive fermions in lattice gauge theory, 
Phys. Rev. D 55 (1997) 3933 [arXiv:hep-lat/9604004]. 

[4] M. B. Oktay and A. S. Kronfeld, New lattice action for heavy quarks, Phys. Rev. D 78 (2008) 014504 
[arXiv:0803.0523 [hep-lat]]. 

[5] A. S. Kronfeld, Application of heavy-quark effective theory to lattice QCD I: power corrections, 
Phys. Rev. D 62 (2000) 014505 [arXiv:hep-lat/0002008]. 

[6] J. Harada, S. Hashimoto, K. I. Ishikawa, A. S. Kronfeld, T. Onogi, and N. Yamada, Application of 
heavy-quark effective theory to lattice QCD II: radiative corrections to heavy-light currents, Phys. 
Rev. D 65 (2002) 094513 [arXiv:hep-lat/01 12044]; Erratum ibid. 71 (2005) 019903. 

[7] J. Harada, S. Hashimoto, A. S. Kronfeld, and T. Onogi, Application of heavy-quark effective theory to 
lattice QCD III: radiative corrections to heavy-heavy currents, Phys. Rev. D 65 (2002) 094514 
[arXiv:hep-lat/01 12045]. 

[8] H. W. Lin and N. Christ, Non-perturbatively determined relativistic heavy quark action, Phys. Rev. D 
76 (2007) 074506 [arXiv:hep-lat/0608005]. 

[9] J. A. Bailey et al. [Fermilab Lattice and MILC Collaborations], The B — > n(v semileptonic form 
factor from three-flavor lattice QCD: a model-independent determination of \V u b\, Phys. Rev. D 79 
(2009) 054507 [arXiv:08 11.3640 [hep-lat]]. 

[10] R. T. Evans, A. X. El-Khadra, and M. Di Pierro [Fermilab Lattice and MILC Collaborations], A study 
of the B s -B s mass and width difference in 2+1 flavor lattice QCD, in proceedings of Lattice 2006, 



PoS (LAT2006) 081 



[11] C. Bernard et al. [Fermilab Lattice and MILC Collaborations], B and D Meson Decay Constants, 



in proceedings of Lattice 2008, PoS (LATTICE 2008)278 [arXiv:0904.1895 [hep-lat]] 



[12] C. W. Bernard et al. [MILC Collaboration], The QCD spectrum with three quark flavors, Phys. Rev. D 
64 (2001) 054506 [arXiv:hep-lat/0 104002]. 

[13] C. Aubin et al. [MILC Collaboration], Light hadrons with improved staggered quarks: approaching 
the continuum limit, Phys. Rev. D 70 (2004) 094505 [arXiv:hep-lat/0402030]. 

[14] C. Bernard, J. Laiho, and R. S. Van de Water, private communication. The most salient formulae can 
be found in Ref. 



[15] R. T. Evans, A. X. El-Khadra, and E. Gamiz [Fermilab Lattice and MILC Collaborations], 

A determination of the B® and B® mixing matrix elements in 2+1 lattice QCD, in proceedings of 



Lattice 2008, PoS (LATTICE 2008) 052 



[16] R. T. Evans, E. Gamiz, A. X. El-Khadra, and M. Di Pierro [Fermilab Lattice and MILC 
Collaborations], A determination of the B® and B®, mixing parameters in 2+1 lattice QCD, 



in proceedings of Lattice 2007, PoS (LATTICE 2007)354 [arXiv:0710.2880 [hep-lat]] 



7 



